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Summary. The precision reached by the recent CMB measurements gives new 
insights into the shape of the primordial power spectra of the cosmological pertur- 
bations. In the context of inflationary cosmology, this implies that the CMB data are 
now sensitive to the form of the inflaton potential. Most of the current approaches 
devoted to the derivation of the inflationary primordial power spectra, or to the in- 
flaton potential reconstruction problem, rely on approximate analytical treatments 
that may break down for exotic models. In this article, we numerically solve the 
inflationary evolution of both the background and all the perturbed quantities to 
extract the primordial power spectra exactly. Such a method solely relies on General 
Relativity and linear perturbation theory. More than providing a tool to test ana- 
lytical approximations, one may consider, without complications, the treatment of 
non-standard inflationary models as those involving several fields, eventually non- 
minimally coupled to gravity. 

The usefulness of the exact numerical approach to deal with CMB data is illus- 
trated by analysing the WMAP third year data in the context of single field models. 
For this purpose, we introduce a new inflationary related parameter encoding the 
basic properties of the reheating era. This reheating parameter has significant ob- 
servable effects and provides a self-consistency test of inflationary models. As a 
working example, the marginalised probability distributions of the reheating and 
potential parameters associated with the small field models are presented. 



1 Motivations 

The inflationary paradigm is currently passing all the tests raised by the 
so-called high precision cosmology measurements ili. Although this suggests 
that the existence of a quasi-exponential accelerated era in the early universe 
may be view as a standard lore, one has to keep in mind that almost all 
the inflationary field models lasting more than sixty efolds and leading to an 
almost scale-invariant power spectrum for adiabatic scalar perturbations may 
do the job. It is therefore of both theoretical and observational interest to look 
for inflationary properties that are, or will be in a foreseeable future, significant 
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enough in the data to allow disambiguation between the different models. 
Many works are devoted to this task ranging from the details of the reheating 
era to the search of a theoretical embedding of inflation in supersymmetry 
or string theory [21 El SI E] • In the following, we will be interested in the 
model disambiguation problem through the cosmological perturbations and 
the Cosmic Microwave Background (CMB) anisotropics. 

Among the analytical tools available to study inflation in the cosmological 
context, the so-called slow- roll approximation provides analytical expressions 
for both the field evolution and the primordial scalar and tensor power spectra. 
It relies on an order by order expansion in terms of the so-called Hubble-flow 
functions £i(n), where n = ln(a/<Xini) is the number of efolding from the be- 
ginning of inflation and a the Friedman-Lemaitre- Robertson- Walker (FLRW) 
scale factor [6]. The Hubble flow fonctions arc defined from the Hubble pa- 
rameter H(n) by 

dlnff dlne^ 

£i = 1 , = —\ ■ (1 

an an 

If the underlying field model is such that these functions remain small at the 
time where the length scales of cosmological interest today leave the Hubble 
radius, then the scalar and tensor power spectra can be Taylor expanded 
around a given pivot wavenumber k„. At first order, one gets [71151 [51 [TU] 



V c (k) 



for the scalar modes and 

2n 2 H 2 



1 - 2 (C + 1) e 1 - Ce 2 - (2e 1 - e 2 ) In * 



V h {k) = 



l-2(C+l)e 1 -2e 1 ln^- 



(2) 



(3) 



for the tensor modes. In Eqs. ([2]) and ([3]), k 2 = Sir /mi is the gravitational 
coupling constant and C is a constant coming from the Taylor expansion 
(C ~ —0.73). The Hubble parameter and the two first Hubble flow functions 
are evaluated at N* = n cn d — n*: the number of efold before the end of inflation 
at which the pivot length scale crosses the Hubble radius: fc* = a(N^)H(N^). 
From these power spectra, assuming the conservation of the comoving cur- 
vature perturbation after Hubble exit (k < aH), the CMB anisotropies in- 
duced by the scalar and tensor perturbations can be derived and compared 
with the data. Using Markov-Chains-Monte-Carlo (MCMC) methods, one 
can extract constraints on the power spectra parameters, namely e±, e 2 and 
P* = k 2 H 2 1 '(8tt 2 : €i) . Assuming a flat FLRW universe, the WMAP third year 
data lead to the posterior probability distributions plotted in Fig. [T] 

A great advantage of the slow-roll approach is that one does not need to 
specify an explicit model of inflation jTOj [T2l [13l [M] . The constraints verified 
by e 1 and e 2 apply to all (single field) inflationary models verifying the slow- 
rolling conditions e x <C 1 and e 2 -C 1 (see however Refs. [HI [16]). However, 
using these results for a given model of inflation requires the knowledge of 
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Fig. 1. WMAP third year data constraints on the first order Hubble flow parameters 
t\(N*) and e 2 (N t ) [11]. The two dashed contours are the 68% and 95% confidence 
intervals associated with the two-dimensional marginalised posterior probability dis- 
tribution. The solid curve corresponds to a scale invariant power spectrum whereas 
the short segments are the slow-roll predictions for the large field models V(tp) oc ip p . 
Note that the model predictions are not "dots" in the plane (ei, e 2 ) due to their de- 
pendence with respect to N* ■ Indeed, due to uncertainties on the reheating era, the 
efold N* for which the observable pivot scale leaves the Hubble radius during in- 
flation is not known (see Fig. [5J). However, under reasonable assumptions, one may 
assume 40 < N t < 60 thereby leading to a "segment" in the plane (e 1 ,e 2 ). 



to determine the associated theoretical values of ej(A^). As can be seen 
in Fig. the value of iV* depends on the number of efold iV re h during which 
the universe reheated before the radiation era. The reheating era depends 
on the microphysics associated with the decay of the inflaton field whose 
complexity renders the determination of N Teh difficult [T71HS1 [TUl 1201 12H 1221 
[23]. However, under reasonable assumptions, it has been shown in Ref. [25] 
that typically 40 < < 60, although these bounds may vary by a factor 
of two for extreme models. For the large field models represented in Fig. [T] 
the uncertainties on the reheating blur the theoretically predicted values of 
the slow-roll parameters (see the short segments in Fig. [T]). Notice that the 
resulting errors in the remain small compared to the current CMB data 
accuracy, but this is not necessarily the case for other models, as for instance 
the small field models discussed in the following. The problem is expected to 
become even more significant with the next generation of more accurate CMB 
measurements. 
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Fig. 2. Sketch of length scale evolution in inflationary cosmology. The horizontal 
axis is the number of efold while the vertical axis represents a logarithmic measure 
of lengths. The cosmological stretching of the observable wavelengths is represented 
by the three blue dashed lines. The evolution of the Hubble radius is represented by 
a solid line from inflation to the matter era. In between, the reheating era connects 
the end of inflation to the radiation era. Although the redshift of equality is known, 
the uncertainties existing on the reheating lead to uncertainties on the redshift at 
which the observable wavelength today have left the Hubble radius during inflation. 
As a result, even for a given model of inflation, the resulting power spectra can be 
significantly different if the number of efold during reheating is changed. This is 
illustrated on this plot by the observability of an inflaton potential feature that may 
accordingly be observable or not. 



Another difficulties that may show up with the slow-roll approach concern 
the existence of features in the inflaton potential. Although the current data 
support the almost scale invariance of the primordial power spectra, the pres- 
ence of sharp localised deviations is still allowed by the data and might even be 
favoured [TH [25l [27l [28l [29J [30J |3T] . In the framework of field inflation, fea- 
tures in the power spectra generically result from transient non-slow-rolling 
evolution associated with sharp features in the inflation potential. In these 
cases, deriving analytical approximations for the perturbations require the use 
of more involved methods [32, 33, 34, [33] . Let us also stress that one of the key 
ingredient rendering analytical methods attractive is the conservation of the 
comoving curvature perturbation £ on super-Hubble scales. As illustrated in 
Fig. [2] this allows one to identify the scalar and tensor power spectra deep in 
the radiation era, which seed the CMB anisotropics and structure formation, 
to the ones derived a few efolds after Hubble exit during inflation. However, as 
soon as inflation is driven by more than one field, the existence of isocurvature 
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modes that may source £ after Hubble exit requires that the modes evolution 
should be traced till the end of reheating [3S] . Although analytical methods 
can still be used, their use is restricted by their domain of validity [5711551155] . 

The previous considerations suggest to use numerical methods to directly 
compute the inflationary perturbations and deduce the primordial power spec- 
tra. Numerical integrations in inflation are not new and have been used to test 
the validity of analytical approximations, or to derive the shape of the power 
spectra for some particular models [40l EH [42J 03J HU [45] . However, with the 
advent of MCMC methods in cosmology, one may be interested in merging a 
full numerical integration of the perturbation during inflation with the CMB 
codes such as CAMB and COSMOMC [46l [47] . The advantage of such an ap- 
proach is that the inflationary parameters become part of the cosmological 
model under scrutiny. When compared with the data, one may expect to get 
consistent marginalised constraints on both the parameters entering the infla- 
ton potential and the usual cosmological parameters describing the radiation 
and matter content of the observed universe. From a Bayesian point of view, 
this is the method that should be used when one is interested in assessing 
the probability of one model to explain the data (statistical evidence). More- 
over, since the underlying approximation is just linear perturbation theory, one 
may consider without additional complications the treatment of non-standard 
models as those involving several fields. Let us mention that our objective is 
to use numerical methods in a well defined theoretical framework. For a given 
model, we use the CMB data to constrain the theoretical parameters and 
discuss the overall model ability to explain the data. Numerical (non-exact) 
methods have also been used in the context of the potential reconstruction 
problems where the goal is to constrain the shape of the inflaton potential 
along the observable window [48j[49l[50] (see Fig. [2]). 

However, as discussed before, to solve the cosmological perturbations from 
their creation as quantum fluctuations during inflation to now, it is necessary 
to model the reheating era. In fact, the importance of reheating for inflation 
is very similar to the importance of the reionisation for the CMB anisotropies. 
Although reionisation of the universe is a complex process, its basic effects on 
the CMB anisotropies can be modelled through the optical depth r. Similarly, 
we will show that the basic effects induced by the reheating on the inflationary 
perturbations may be taken into account through a new parameter In R which 
has not been considered so far. 

The plan is as follow. In a first section, the theoretical setup is introduced. 
We use the sigma-model formalism which allows an easy implementation of 
any scalar field inflation models in General Relativity and multi-scalar tensor 
theories. The equations of motion for the background fields and their pertur- 
bations are presented in the first section. Their numerical integration is the 
subject of the second part. After having introduced our modelisation of the 
reheating era, the last section illustrates the usefulness of the exact numerical 
method by an analysis of the third year WMAP data in the context of the 
small field models 
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2 Multifield inflation 



It is out of the scope of this work to deal with all the inflationary models pro- 
posed so far. However, four-dimensional effective actions associated with many 
inflation models, and especially the ones being embedded in extra-dimensions, 
share the common feature that they involve several scalar fields that may be 
non-minimally coupled to gravity. This is for instance the case for the moduli 
associated with the position of the branes in various string-motivated infla- 
tion models [51] ■ It is therefore convenient to consider an action that may 
generically drive the dynamics of both minimally and non-minimally coupled 
scalar fields, as in the sigma-model [52] [53] [54] [55] . 



2.1 Sigma-model formalism 



Denoting by J 7 ""^ 1 ) the n„ dimcnsionless scalar fields living on a sigma-model 
manifold with metric £ab{^ c ), we consider the action 



S 



1 

2k 2 



R ~ Ug^d^dvF* - 2V{T C 



-g d 4 £, 



(4) 



where g^ is the usual four dimensional metric tensor of determinant g, R the 
Ricci scalar and V the field potential. 

For instance, if <p = J-^ / n and In = 1, this action describes a unique 
minimally coupled scalar field. In this case, the associated potential is U(ip) = 
V{ip)/n 2 and we recover the standard form 



1 



S = — r / RyJ—g d X 



2k 2 



U(<p) 



-gd A x. (5) 



Another example is provided by the models of brane inflation where the in- 
flaton field (p lives on the brane and a bulk field \ m the four-dimensional 
effective action couples to gravity in a non- minimal way [561 157[ 15 lj . With 
n a = 2, tp = J 7 ^ /k and x = 3~( 2 \ Eq. ([4]) can be recast into 



-gd x 



-gd x, 



(6) 



for 



£ ab =diag(A 2 , 1), V(<p, X ) = W( X ) + k 2 A 4 (x)U(<p). 



(7) 

This action describes the dynamics, in the Einstein frame, of the field tp 
evolving in a potential U in a scalar-tensor theory of gravity where x 
is the scalar partner to the graviton [58]. The conformal function A 2 (x) 
and the self-interaction potential W{x) depend on the brane setup consid- 
ered [59] [60] [61] [62]. In the general case, Eq. (J4| can be used to describe 
multifield inflation in a multi-scalar tensor theory of gravity. 
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Differentiating the action (HJ) with respect to the metric leads to the Ein- 
stein equations 

G>„ = S^ Vl (8) 

with the source terms 

= i ab Sf v - g^V, (9) 

where 



S$ = d^T a d v T h - -g, lv d p T a d?T h . (10) 



1 

Similarly, the fields obey the Klein-Gordon-like equation 

UT C + g^T c ah d p T a d v T h = V c , (11) 
where Y denotes the Christoffel symbol on the field-manifold 

?ab = ^ (tda,b + idb,a ~ Lb,d) , (12) 

and V c should be understood as the vector-like partial derivative of the po- 
tential 

BV 

yc ={ cdy d=e cd_ (13) 

2.2 Background evolution 

In a flat (FLRW) universe with metric 

ds 2 = g^dx^dx" = a 2 (r]) (-d V 2 + S lJ dx i dx j ) , (14) 

r\ being the conformal time and i and j referring to the spatial coordinates, 
the equations of motion ([5]) and (fTTj) simplify to 

m 2 = ^labT' 1 '^' + a 2 V, (15) 

2H' + H 2 = -\^abT a 'T b ' + a 2 V, (16) 

T c " + T c ah T al T hl + 2HT CI = -a 2 V c , (17) 

where a prime denotes differentiation with respect to the conformal time and 
H = aH is the conformal Hubble parameter. In terms of the efold time variable 
n, the field equations can be decoupled from the metric evolution and one gets 

V 

H 2 = j— , (18) 

3 -r 2 v 
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a dot being a differentiation with respect to n. We have introduced a velocity 
field 



a = \JlabT a T h . (21) 

In fact, a is the so-called adiabatic field introduced in Ref. [36] which describes 
the collective evolution of all the fields along the classical trajectory. From 
Eqs. (fl9|) and Eq. (|20|) one may determine its equation of motion 

a" + 2Ha' + a 2 V a = 0, (22) 

with V a = u°V c and where the u a are unit vectors along the field trajectory: 

u« = i— = (23) 
a g 

In the Einstein frame, inflation occurs for d 2 a/dt 2 > 0, or in terms of the first 
Hubble flow function, foiQ 

* — § - |*« < 1 < 24 > 

The multifield system induces an accelerated expansion of the universe if 
the resulting adiabatic field velocity a remains less than y/2. According to 
Eq. (|20[). the term in 1/(3 — <r 2 /2) acts as a relativistic-like inertia for the fields 
evolution and thus & < y6 (for a positive potential) . In this equation, the first 
term on the left hand side may be interpreted as a covariant acceleration on 
the curved field manifold, the second as a constant friction force and the 
right hand side as a driving force deriving from the potential lny. In fact, 
we recover the well-known attractor behaviour of the inflationary evolution: 
whatever the initial fields velocity, the friction term ensures that the terminal 
velocity of the fields will be, after a transient regime, 

Analytical integration of the previous expression is at the basis of the slow- roll 
approximation when the effective potential In V is flat enough. In the general 
case, the driving force is always pushing the fields towards the minimum of 
In V. Let us note that this is why the monomial potentials V oc tp p are actually 
"flat" for the large field values: d In V/d(p = p/ip which goes to zero for <p — > 00. 



2.3 Linear perturbations 
Scalar modes 

In the longitudinal gauge, the scalar perturbations (with respect to the rota- 
tions of the three-dimensional space) of the FLRW metric can be expressed 
as 

1 Notice that this does not imply that the universe is accelerating in the string 
frame and one has to verify that there are enough efolds of inflation to solve the 
homogeneity and flatness issues in that frame [551 1631 164] . 
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ds 2 = a 2 [- (1 + 2#) d?] 2 + (1 - W) ll3 Ax l Ax j ] , (26) 

where $ and W are the Bardeen potentials. With 8T a the field perturbations, 
the Einstein equations perturbed at first order read 

m&' + (H' + 2U 2 ) <P - AV = -\ 2 i ah T al &T hl 

' f \tab,cT al T b ' + a 2 v)j5T c , (27) 
&' + H& = \^ ah T a '8T h . s (28) 

v" + zm>' + in 1 + 2H 2 ) & = -iab^'ST 1 " 

where use has been made of <P = from the perturbed Einstein equation with 
i 7^ j. For flat spacelike hypersurfaces, 

A = 5 l3 d t d v (30) 

Similarly, the perturbed Klein-Gordon equations read 

5T C " + 2T c ah T al ST hl + 2H8T cl 

+ (r c ab,dF al F h ' + a 2 vi - t a i ab4 a 2 V b ^j 8T d 

- A8T C = W'T cl - 2&a 2 V c . (31) 

As discussed in the introduction, if there is more than one scalar field involved, 
the entropy perturbation modes can source the adiabatic mode even after 
Hubble exit. The equation governing the evolution of the comoving curvature 
perturbation £ can be obtained from Eqs. (|27p to (f2"9"|) , using the background 
equations. Firstly, the Bardeen potential verifies 

+ ' + (2H 1 + AH 2 ) & - A& = -a 2 V c 8T c . (32) 

Using the geometrical definition for the comoving curvature perturbation |13j 

t = n U _ n2 (& + H$), (33) 

Eq. (US]) yields 

( = * + H 6 -^, (34) 
a 

where the adiabatic perturbation 8a is also the resulting perturbation of all 
fields projected onto the classical trajectory [see Eqs. |2T|) and (|23l) ]: 
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5a = ab = u a 5T a . (35) 

a' 

The dynamical equation ([3^| now exhibits couplings between the adiabatic 
and entropy modes 

C -™*-™(jVjr-^±™*), (36) 
a' A u \ a a J 

which can be recast into 

<?=™t®-™L\*VJT*. (37) 
g a ■ L 

The orthogonal projector is defined by 

J-a&= lab ~ Vab, (38) 

where r]^ = u a Uf, is the first fundamental form of the one-dimensional man- 
ifold defined by the classical trajectory [65] . Clearly, the comoving curvature 
perturbation on super-Hubble scales for which A'F ~ is only sourced by 
the entropy perturbations defined as the projections of all field perturbations 
on the field-manifold subspace orthogonal to the classical trajectory. If, on 
the contrary, there is a single field involved during inflation, then these terms 
vanish and we recover that £ remains constant after Hubble exit. 



Tensor modes 

In the Einstein frame, the scalar and tensor degrees of freedom are decoupled. 
Therefore, the equation of evolution for the tensor modes remains the same 
as in General Relativity. For a flat perturbed FLRW metric 

ds 2 = -a 2 drf + a 2 {6 tJ + h tJ ) dx l dx 3 , (39) 

where hij is a traceless and divergenceless tensor 

.Y'l,,, d-'ihh,, 0. (40) 

one gets [3 [66] 

h'-j + 2Hh' l0 - Ah l3 = 0. (41) 
Primordial power spectra 

The initial conditions for the cosmological perturbations require the knowl- 
edge of the two-point correlation functions for all of the observable scalar and 
tensor modes deep in the radiation era. In Fourier space, these are just the 
power spectra associated with the values taken by the adiabatic and entropy 
perturbations at the end of the reheating, i.e. 



The exact numerical treatment ol inflationary models 



11 



V ab = ^[v*{k)Y[v\k)], (42) 

where v a stands for C, or the entropy modes. Similarly, taking into account 
the polarisation degrees of freedom, the tensor power spectrum reads 

V h {k) = \\h{k)\\ (43) 

In the following, we summarise the numerical method used to solve the 
full set of Einstein and Klein-Gordon equations derived in this section. The 
power spectra can then be deduced from Eqs. (|42l) and (|43l) by pushing the 
integration till the end of the reheating. 



3 Numerical method 

3.1 Integrating the background 

As suggested by the form of Eqs. (fl8| to (|20|) . it is convenient to use the 
number of efold n as the integration variable. In fact, the background evolution 
only requires the integration of the fields equation of motion (|20p . From the 
Cauchy theorem, the solution is unique provided all the T a {Q) and J- a (0) 
are given at n — 0. Plugging the solutions for F a (n) into Eq. (fT8|) uniquely 
determine the Hubble parameter and thus the geometry during inflation. 



Initial conditions 



However, as previously mentioned, the attractor behaviour induced by the 
friction term erases any effect associated with the initial field velocities after 
a few efolds. This is the very reason why initial conditions in inflation are 
essentially related to the initial field values only. On the numerical side, the 
attractor ensures the stability of almost all forward numerical integration 
schemef0. In the following, we have used a Runge-Kutta integration method 
of order five and the initial field velocities have been chosen on the attractor 
by setting 



(44) 



The robustness of the attractor during inflation may be quantified by com- 
paring the numerical solutions obtained from various arbitrary choices of the 
initial field velocities, at fixed value of !F a (0). As an illustration, we have plot- 
ted in Fig. [3] the efold evolution of .F a (n) in a brane inflation model involving 
three scalar fields, two of them are non-minimally coupled to gravity and rep- 
resent the position of two branes in a five-dimensional bulk (see Ref. [SS]). As 
can be seen on this plot, all the fields are on the attractor after a few efolds. 



2 as well as the instability of backward integrations. 
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Fig. 3. Evolution of the field velocities T a = {%, 0, ip} in the boundary inflation 
model of Ref. |55| which involves three fields, two of them being non-minimally 
coupled to gravity. The dotted curves are the solutions obtained by setting the 
initial field velocities on the attractor according to Eq. (|44p whereas the solid and 
dashed curves are the solutions obtained from a random choice of the initial field 
velocities (ensuring however H 2 > 0). The first Hubble flow parameter in the string 
frame is plotted in the lower panel. Since this is also the adiabatic field velocity 
squared, the acceleration properties of the universe after a few efolds do not longer 
depend on the initial field velocities [see Eq. (1241) ]. 



End of inflation 



From the above initial conditions the fields evolve toward the minimum of 
the potential InV while the expansion of the universe accelerates as long 
as e x < 1. It would therefore be natural to define the end of inflation by 
the efold n en( j at which e^riend) = 1- However, this is usually not the end 
of the fields evolution since they have not yet reached the minimum of the 
potential. On the contrary, £i(n cn d) = 1 just signals that the kinetic terms 
in Eq. (fH) start to dominate over the potential. Since the expansion factor is 
decelerating for ei(n cnc i) > 1, this late stage evolution takes place during a 
few efolds and the fields rapidly reach the minimum of the potential. In the 
standard picture, the fields oscillate around the minimum of the potential and 
decay through parametric resonances into the relativistic fluids present during 
the radiation era [TH1 [57] (see Fig. 0]) The details of the reheating process 
are very model dependent and require the knowledge of all the couplings 
between the inflaton and the standard model particles [2]. This implies that 
the number of efolds the universe reheated also depends on the model at hand. 
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Another complications come in multifield inflationary models in which the end 
of inflation and the reheating may be triggered by tachyonic instabilities. In 
these cases the condition e 1 (n cnc i) = 1 is not longer relevant and one should 
rather introduce a limiting field value .Fend to classically define the beginning 
of the reheating era (and the end of inflation) . 

Following the previous discussion, our phenomenological approach to the 
end of inflation is to assume an instantaneous transition to the reheating era. 
The efold n en( j at which the transition occurs is either determined by the 
condition ei(n en d) = 1 or when the relevant field cross a limiting value Fend! 
the choice being made according to the inflation model we are interested in. 
The value of F cn d may be given by the underlying microphysics or considered 
as an additional parameter of the inflation model. For instance, as can be seen 
in Fig. HJ it is convenient to define the end of inflation for the large field models 
by e l = 1: the field evolution afterwards, oscillations around the minimum of 
the potential and subsequent decay, is supposed to be part of the reheating 
stage. 

Knowing the fields value T a {n) from n = to n = n en d; the background 
geometry is given by Eqs. (|18|) and (|19p and we can now numerically inte- 
grate the linear perturbations on the same efold range. As discussed in the 
introduction, the link with the cosmological perturbations observed today still 
requires a reheating model that will be introduced in Sect. 14.11 
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3.2 Integrating the perturbations 

As for the background, we have chosen to integrate the linear perturbations 
in efold time. Focusing on the scalar perturbations, their dynamics is driven 
by the Einstein and Klein-Gordon equations given in Sect. 12.31 These ones are 
however redundant due to the stress energy conservation already included in 
the Bianchi identities. As a result, it is necessary to integrate only a subset 
of Eqs. (|27p to (|3"Tj) . Although the Bardeen potential & could be explicitly 
expressed in terms of the field perturbations 8T a only, such an expression 
is singular in the limit k — > and e ± — > 0, which is not appropriate for a 
numerical integration (see below). It is more convenient to simultaneously 
integrate the second order equations (|3"Tj) and (|3"2")) . Recast in efold time, they 
read 

8T- + (3 - ejSfc + 2T c ab T-8T b + (r^ lbd T-T b + ^ a U,d^j 

+ ^Jp S ^ = 4 ^ C " 2 ^ ( 45 ) 

* + (7 - *)* + (a^r + jjp) * = -^ C - ( 46 ) 

The constraint equations (|2"T|) and (|2"5)l being first integrals of the above equa- 
tions, there is still an integration constant that should be set to restore the 
equivalence to the full set of Einstein and Klein-Gordon equations. This one 
can be fixed by choosing the appropriate initial conditions at n = n- lc for the 
Bardeen potential \P. Setting all the 8T a (ni c ) and 8j- a (n lc ), the initial con- 
ditions for the Bardeen potential are indeed uniquely given by Eq. (f2"T|) and 
([2"8jl . In efold time, one gets 



a 2 H< 



$ = \labT a 8T h - V , 

(47) 

these expressions being evaluated at the initial efold time. As a result, the 
linear perturbations of both the fields and metric are uniquely determined by 
the initial conditions 5T a {n- lc ) and 8T a (n lc ). As discussed in the next section, 
the initial conditions are set on sub-Hubble scales for which k ^> aH ensuring 
the regularity of Eqs. (j47|) . 



Quantum initial conditions 



In the context of single field inflation, the initial conditions for the linear per- 
turbations are given by the quantum fluctuations of the field-metric system on 
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sub-Hubble scales k — > oo. In this limit, the perturbations decouple from the 
expansion of the universe and a field quantisation can be performed along the 
lines described in Refs. [EHH]. The canonically normalised quantum degrees 
of freedom are encoded in the Mukhanov-Sasaki variable 

Q = 5a+—& = 5a+ y/&, (48) 
n 

with a — = Kip for single field models. In terms of Q, the equation of 
motion ([3"Tj) can be recast into 

(aQ)" + 

showing that in the small scales limit k —> oo, the quantity aQ follows the 
dynamics of a free scalar field. Assuming a Bunch-Davies vacuum, aQ has a 
positive energy plane wave behaviour on small scales and in Fourier space one 
gets 

e -ikr/ 

lim aQ(ri) — k — -j^. (50) 

This solution uniquely determines the subsequent evolution of the perturba- 
tions during inflation and will be our starting point for the numerical inte- 
gration. According to Eq. (|4"T|) , provided the initial conditions are set in the 
limit k/H — > oo, Eq. (fSTJf is also the small scales behaviour of the rescaled 
adiabatic field perturbations aSa. 

For a multifield system the previous results can be generalised and the 
quantum modes identify with the adiabatic perturbations 5a together with 
the canonically normalised entropy modes introduced in Sect. 12.31 In fact, if 
the original fields T a are already canonically normalised, i.e. i a \, = S a b, arid 
independent dynamical variables in the small scales limit, the adiabatic and 
entropy perturbations can be obtained from the original field perturbations by 
local rotations on the n CT -dimensional field manifold 36, 38, 55 . In this case, 
denoting by 5^ a the adiabatic and entropy perturbations, with the convention 
£<jW = 5a 7 one has 

6s a = M a b {T c )8T b , (51) 

where M. is an instantaneous rotation matrix, = M. T = M^ 1 , depending 
on the background quantities only. Under these assumptions, the quantum 
modes are independent in the small scales limit and their two point correlators 
reduce to 

(<k a *(fc)<k b (fc')) k = n S ab V,{k)S(k - k'), (52) 

being the free field power spectrum given by the square modulus of Eq. (|50p . 
Consequently, all the correlators between the original field perturbations in- 
herit these initial conditions: 



aQ = 0, 



(49) 
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(5T a *{k)5T\k')) = (M -1 )"* (M~ 



l ) ,(d^ c *Sq d ) = 8 ab V,(k) S(k - k') 



(53) 

The previous results can be generalised for the sigma-models with a diagonal 
metric £ a b non equal to the identity by the transformation 8T a — > ^/£ a aSJ- a 
(no summation). In the general case, the transformation matrix .M would mix 
all the fields and the cross correlators. However, since it is always possible to 
diagonalise £ a b through a field redefinition, we will now assume without lost 
of generality that £ a t is diagonal. 

Defining the normalised quantum modes by 



£ = aV2eV 2 k 3 / 2 5f a , 



(54) 



with l aa evaluated along the background solution, in the small scales limit 
Eq. ([50]) can be recast into the initial conditions 



Ms' 

/i q — Kk 7 - = 

fc>W k k»H 



(55) 



up to a phase factor. In terms of the field perturbations, the initial conditions 
in cfold time therefore read 



k 3/2 ST a 
k 3/2 5T a 



ic a a ond v / £ aa (n ic ) ' 
nk a n e Tlcnd ~" ic 



1o Clend y/£ aa (rii c ) 



, 1 ^aa(^ic) i ^ 
2 ^aa (^ic) 



(56) 



The efold n lc at which these initial conditions should be set has not been 
specified yet. In fact, the limit k/TL — >■ oo would correspond to the infinite 
past and does not make sense for non-eternal field inflation model^]. However, 
by definition of n = In (a/a- m i), the condition k/a 3> H is already satisfied a 
few efolds before Hubble exit. For all inflation models lasting more than iV* 
efolds, it would be natural to set the initial conditions for the perturbations 
at the beginning of inflation n ic = (see Fig. [5]). 

Choosing n ic = is however not appropriate for a numerical integration. 
Indeed, according to the initial values of the background fields, the total num- 
ber of efolds n cn d can be much greater than N*. In such cases, most of the 
computing time for the perturbations would be spent into the deep sub-Hubble 
regime for which the modes behave as free plane waves. It is rather more con- 
venient to set the initial conditions "closer" to the time nk at which a given 
mode cross the Hubble radius k — H(nk). Following Ref. [3D], a simple choice 
is to define nj C for each mode according to 



3 Eternal inflation may occur when the quantum fluctuations on Hubble length 
scales become dominant over the classical field evolution and the semi-classical 
approach used here would no longer be valid, at least in the self-reproducing 
regime. 
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k 



= a 



(57) 



H(n ic ) 



C q being a constant verifying C q 3> 1 and characterising the decoupling limit. 
Strictly speaking, this choice introduces small trans-Planckian-like interfer- 
ences between the modeifl which remain however negligible provided C q is big 
enough [68] . 

Mode integration 

For each perturbation mode of wavenumber k, Eqs. (|45p and (I46[) are numer- 
ically solved by setting the initial conditions at the efold n ic , solution of 
Eq. ([5T)l . In order to significantly speed-up the numerical integration, instead 
of using the already computed background solution it is more convenient to 
integrate both the background and the perturbations simultaneously. This 
can be done along the following steps. 

Firstly, Eq. ([57]) is solved to determine ni c (fci) for the largest wavelength 
mode k\ we are interested in. The background equations (|20p and are then 
integrated from n = to n = n ic (fci). At that efold, Eqs. (HOI), (@5J and gBJ) 
are simultaneously integrated till the end of inflation at n = n on( j. This pro- 
cess is iterated for each of the ki > ki_\ mode wanted. However, to speed-up 
the integration, it is enough to re-integrate the background from n\ c (ki-\) to 
ri\ c (ki) rather than from n — before switching on the perturbations (see 
Fig. [5]). Such an integration gives the value of all the field and metric pertur- 
bations at the end of inflation n = n enc i. In principle, the power spectra can 
then be deduced by using Eq. (|42| . 

Primordial power spectra 

For a multifield system, since the perturbation modes are supposed to be 
independent deep under the Hubble radius, they can be considered, from a 
classical point of view, as independent stochastic variables. As a result, the 
power spectra at the end of inflation are no longer given by Eq. ([42]) but 
should be computed as [13] 



4 In the free field limit, a more accurate choice for ni c is kr](n\ c ) = C q . This def- 
inition maintains the phase factor in Eq. (|50[) independent of k on the initial 
hypersurface. 

For direct numerical integrations, each step requires various forward and back- 
ward evaluations of the background functions. If these ones are not analytically 
known but precomputed, one has to use spline and interpolation methods which 
are heavily time-consuming. 




(58) 



m— 1 
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Fig. 5. Sketch of the numerical integration for the scalar perturbation in the m 2 ip 2 
single field model. In the top frame, the initial conditions for each mode are set at 
the efold n(fc;) solution of Eq. (|57fl . The perturbations (and the background) are 
then integrated till the end of inflation at n = n en d- The bottom frame represents 
the efold evolution of the Mukhavov-Sasaki variable Q for the two corresponding 
modes [see Eq. (|48[) ]. Notice that only £ = Q/\/2e 1 is conserved after Hubble exit. 



where, as before, v a stands for the observable perturbations one is interested 
in. The v a can be the field perturbations themselves but it is more custom- 
ary for CMB analysis to use the comoving curvature perturbation £ and the 
rescaled entropic perturbations 8z, a ja (a > 1). The index "to" in Eq. 
refers to the n a independent initial conditions obtained by setting only one 
perturbation mode /i™ in the Bunch-Davies vacuum at n; c , the other \x^ m 
vanishing. Notice that we have not explicitly written the entropy modes since 
various definitions are used in the literature. The definition of the entropy 
modes through the standard orthogonalisation procedure along the field tra- 
jectory can be found in Refs. [351 [35] an d has the advantage to give canonically 
normalised perturbations. Another definitions introduce a reference field and 
define the entropy perturbations to be the relative perturbations of the other 
fields with respect to it. These differences come from the fact that one has to 
specify how the fields decay after inflation to know between which cosmolog- 
ical fluids entropy perturbations may exist. For instance, if all the cosmolog- 
ical fluids observed today are produced by the decay of one field only, then, 
although entropy perturbations exist during inflation, they are usually not 
observable afterwards (69} [70] ■ Their only effect would be to break the conser- 
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vation of £ on super-Hubble scale thereby requiring the integration of all the 
perturbations till the end of inflation to determine the £ power spectrum. 




Fig. 6. Power spectra of the comoving curvature perturbation £ at the end of 
inflation from the first order slow-roll approximation (dotted line) and an exact 
numerical integration (V tx ip 2 ). The wavenumbers are expressed through n^, the 
efold at which the mode k crosses the Hubble radius: k = TL(nk) (see Fig. [5}. 



As an illustration, Fig.[S]shows the exact numerical power spectrum for the 
comoving curvature perturbation £ for the single field chaotic model V oc ip 2 . 
As can be seen on this plot, although the exact power spectrum differs from 
its first slow-roll approximated version given by Eq. the differences remain 
small on a 10 efold observable range. Another example involving entropy per- 
turbations is plotted in Fig. [7] for a two fields model of inflation. The presence 
of one entropy mode breaks the conservation of £ on super-Hubble scale and 
the so-called consistency check of inflation Vn = 16e x V( [see Eqs. © and 
at zero order]. 

Tensor perturbations 

Since in the Einstein frame the tensor and scalar degrees of freedom are de- 
coupled, the numerical integration of the tensor modes does not present any 
difficulties. The equation of motion (|TT|) can be recast into the equation of 
a parametric oscillator by defining the canonical mode function fi T = k 3 ^ 2 ah 
satisfying the deep sub-Hubble initial conditions ([50)1 for a Bunch-Davies vac- 
uum [8]. The power spectrum Vh is readily obtained by evaluating Eq. (|43|) 
at the end of inflation. 
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n k" "end 

Fig. 7. Violation of the "consistency check" of inflation in a non-minimally coupled 
two-fields model. The action is given by Eq. Q with W = 0, U cc <p 2 and A 2 — 
exp(— ax)- Even for a small value of a — 1/30, the ratio between the scalar and 
tensor power spectra is no longer equal to 16ei. If the two fields were uncoupled on 
all scales, then the power spectrum of entropy perturbations would be the one of a 
free test scalar field Vsi = 'Pf l /16. As can be seen on the plot, this last condition is 
also violated. 



Physical wavenumbers 

Up to now, one may have noticed that all the power spectra have been plotted 
with respect to rif. and not with respect to the values of k. For astrophysical 
purpose, one needs to know the correspondence between the comoving k ap- 
pearing in the above equations and the physical wavenumbers measured today 
fc/oo, whose typical unit is the Mpc" 1 . As it appears in the initial conditions 
(|56p . rendering fc/ao explicit requires the knowledge of ao/a en d, i-e. the red- 
shift z cn d associated with the end of inflation. As discussed in Sect. [1] this 
can only being done if one knows the number of efolds during which the uni- 
verse reheated. Let us also notice that the physics involved in the quantum 
generation of cosmological perturbations appears through these very numbers: 
nk/ao are the wavenumbers measured today, usually of Mpc -1 size, expressed 
in unit of the Planck mass with k = \/8n/m pl [see Eq. (|56|)]. 

4 Application to CMB data analysis 

Fig. [5]makes clear that from the integration of the perturbations described in 
the previous section, their power spectra are known at the end of inflation. 
These primordial correlations then evolve through the reheating, radiation 
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and matter era to shape the universe into its current state. The theory of 
cosmological perturbations precisely predicts how such linear perturbations 
evolve in a FLRW universe from the tightly coupled regime deep inside the 
radiation era to today. As a result, we still have to know how the power spec- 
tra are modified through the reheating era. As already mentioned, reheating 
is very model dependant and a detailed analysis is out of the scope of our 
current approach. Instead, remembering that the objective is to use the CMB 
anisotropics measurements as a probe to get information on the primordial 
correlations, we introduce a basic model of reheating described by some phe- 
nomcnological parameters. 

4.1 Reheating 

Assuming that perturbations on super-Hubble scales are not significantly 
modified till the beginning of the radiation ers^l, the reheating may influence 
the observed power spectra through its effects on z on( j (see Fig. [2]). For in- 
stantaneous transitions between inflation, the reheating era and the radiation 
era [24], one has 

m ^E! = i n ^ +m ^££li + i n ^a ; (59) 

ftp G-rch &eq Q 

where a on d, flroh and a eq are respectively the scale factor at the end of inflation, 
at the end of reheating and at equality between the energy density of radiation 
and the energy density of matter. The redshift of equality can be expressed 
in terms of the density parameter of radiation today f2 r &d and the Hubble 
parameter today H . Moreover, during the radiation era p oc a , and Eq. (|59|) 
can be recast into 

In = hi - 7 In (^V-ch) + ~ hi ( \/ZQ X3A kH^] , (60) 
a Q a rch 4 '2V / 

where p re h denotes the total energy density at the end of the reheating era. 
It is clear that the first two terms depend on the physics involved during the 
reheating. For instance, they would only depend on the energy density at the 
end of inflation for a radiation-like reheating era. This suggests to introduce 
a phenomcnological parameter [llj 

lni2 rad = ln^l--mf^V (61) 

arch 4 \ p cn d / 

From this parameter, the quantity k/TL entering the equations of motion for 
the perturbations (I45[) and (|46|) can be evaluated in terms of the fc/eto values 
measured today 

6 Although such an assumption is motivated by the fact that the physical processes 
involved during reheating are sub-Hubble, this assumption may not longer be true 
in presence of entropy modes. 
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aH a K H(n) (3n iad ) 1/A ' 
and similarly for the initial conditions in Eq. (1561) . 

4.2 CMB anisotropies 

For a given value of i? ra d, the primordial power spectra deep in the radiation 
era are now uniquely determined by the numerical integration described in 
Sect. [3] and can be used as initial conditions for the subsequent evolution of 
the perturbations. The integration of the cosmological perturbations through 
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Fig. 8. Angular temperature power spectrum of the CMB anisotropies in a ylCDM 
universe born under chaotic inflation. Fiducial values of the parameters have been 
used: i? rad = 1, p = 2, kM = 2 x 1(T 3 , fl h h 2 = 0.022, J? c = 0.12, h = 0.7, z IC = 12, 
where h is the reduced Hubble parameter and z le the redshift of reionisation. The 
WMAP third year measurements are represented as blue squares. 



the radiation and matter era, as well as the resulting CMB anisotropies, have 
been performed by using a modified version of the CAMB code |46j . The model 
parameters involved are both the inflation parameters and the usual cosmo- 
logical parameters describing the FLRW model at late time. For instance, for 
a ylCDM universe experiencing large field inflation in its earliest times, there 
are two parameters fixing the potential V(ip) = M 4 ip p , one parameter describ- 
ing the reheating era -R ra d, plus the four cosmological base parameters: the 
number density of baryons i?b, of cold dark matter J? c , the Hubble parameter 
today H a and the redshift of reionisation of the universe z Te [47] . Let us recap 
that we have defined the end of large field inflation by e 1 (</? en( j) = 1- Combined 



The exact numerical treatment of inflationary models 23 

with the existence of the attractor during inflation, this ensures that (/3 en< j is 
fixed by the potential parameters [see Eq. (|55|l]. The resulting angular power 
spectrum for the CMB temperature fluctuations is represented in Fig. [5] for a 
fiducial set of the parameters. 

The next step is to use CMB measurements to constrain the models. For 
this purpose, the parameter space can be sampled by using Markov Chain 
Monte Carlo (MCMC) methods as implemented in the COSMOMC code [47] 
to extract the probability distributions satisfied by the model parameters. In 
the next section, we illustrate such a procedure for the vlCDM model born 
under small field inflation by using the WMAP third year data [711 [72] [73] . 



4.3 WMAP3 constraints on small field models 



Small field inflation can be described by the action §5§ when the potential 
reads 



U(tp) = M 4 



1-1 2 

A* 



(63) 



The inflation model parameters are the energy scale M, the power p and 
the vacuum expectation value scale fi. As can be seen in Fig. [5] inflation 
proceeds for small initial field values and stops when ei(<^ cn d) = 1- Notice 
that the potential in Eq. (|63[) is negative for (f > ji which means that the 
above description is no longer correct. This is however not an issue since it 
occurs well after the end of inflation and the basic effects of the reheating are 
already encoded in our extra parameter i? ra d. 




Fig. 9. The small field potential In V on the left and the first Hubble flow function 
£i on the right. Inflation occurs for £i < 1, for small values of the field. The three 
curves correspond respectively to p = 2, p = 4 and p — 6 from the left to the right 
(from Ref. [II]). 
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Observable parameters 

Since flatness is inherited from inflation, the /1CDM cosmological model is 
described by the density parameters associated with the different cosmolog- 
ical fluids: f?b, fl c plus the Hubble parameter today H . The cosmological 
constant is fixed by Qa = 1 — — However, in order to minimise the 
parameter degeneracies with respect to the CMB angular power spectra, it 
is more convenient to perform the MCMC sampling on the equivalent set of 
parameters f2bh 2 , f2 c h 2 7 the optical depth r and the quantity 9 which mea- 
sures the ratio of the sound horizon to the angular diameter distance (h is the 
reduced Hubble parameter today) [UJ . 

Similarly for the potential parameters, as can be seen in Eqs. ^ and 
([3]), the overall amplitude of the CMB anisotropics is proportional to the 
Hubble parameter squared and thus to the potential V [see Eq. (TT8]) ]. Since 
the amplitude of the cosmological perturbations is a well measured quantity, 
the data may be more efficiently used by directly sampling the primordial 
amplitude of the scalar power spectrum "P* = Vc,(k*) instead of M (fc* being 
a fixed observable wavenumber: fc*/a = 0.05Mpc _1 ). 

However, the numerical method used to integrate the perturbations during 
inflation requires the input of a numerical value for M to predict the value 
of V*. In fact, one can use the trick described in Ref. [IT] : under a rescaling 
V — > sV, the power spectrum scales as V((k) — > s7 3 ^(s 1 / 2 fc) at fixed i? ra d- The 
idea is therefore to integrate the perturbations with an artificial normalisation 
of the potential, for instance M = 1, and then analytically rescale M from 
unity to its physical value that would be associated with V*. The required 
value of s is given by the ratio 'P*/7 : o, M ~ 1) , where r p^ A ~ 1) is the amplitude 
of the scalar power spectrum stemming from the numerical integration with 
M = 1 and evaluated at fc = fc*s -1 ' 2 . Still, it is not really straightforward to 
determine s since both Vq and k change simultaneously. The last subtlety is 
to remark that k = fc„ if instead of considering i? ra a fixed, one considers the 
rescaling of M at fixed In R, with In R defined by 

In R = In i? rad + - In (k Vnd) • (64) 

The parameters R and i? ra d differ only by p cn d, the energy density at the end 
of inflation which is uniquely determined from M, \x and p. It will be therefore 
more convenient to sample the model parameters V* and In R (together with 
/i and p) rather than M and i? ra d- 

Priors 

The prior probability distributions for the base cosmological parameters Q^h 2 , 
f2 c h 2 , t and 9 have been chosen as wide top hat uniform distribution centred 
over their current preferred value [TJ I47j . 
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Concerning the inflaton potential parameters, their priors can be chosen 
according to various theoretical prejudices [TJJ. Since V* is related to the 
energy scale during inflation, a uniform prior has been considered around 
a value compatible with the amplitude of the cosmological perturbations: 
m(10 10, P*) € [2.7, 4.0]. For the scale p associated with the vacuum expectation 
value of if, we have considered two priors. The first includes smaller and larger 
values than the Planck mass: k\i S [1/10, 10]. The second prior is also uniform 
but include values much larger than the Planck mass: up £ [1/10, 100] . Finally, 
a uniform prior is chosen for the power p on [2.4, 10] (p = 2 is a particular 
case, see Ref. [Tip. 

It remains to express our prior knowledge on the reheating parameter In R. 
As previously mentioned, we are assuming that gravity can still be described 
classically which only makes sense if the energy densities involved remain 
smaller than the Planck energy scale, namely for K 4 p nd < 1. On the other 
side of the energy spectrum, the success of big-bang nucleosynthesis (BBN) 
requires that the universe is radiation dominated at that time, thus p rc h > p nuc 
with p nuc ~ 1 MeV 4 (and p enc j > p ro h)- Moreover, we will assume that during 
reheating the expansion of the universe can be described as dominated by a 
cosmological fluid of pressure P and energy density p. In this case, in order 
to satisfy the strong and dominant energy conditions in General Relativity, 
one has —1/3 < P/p < 1 (notice that P/p < —1/3 would be inflation). From 
Eqs. I|6ip and (|64J). the resulting bounds read 

i In («; 4 p nuc ) < In R < -— In (/t 4 /0 nuc ) + jjj In (k 4 /w) , (65) 
and an uniform prior on In R have been chosen in between. 

Results 

The data sets used to constrain the small field ylCDM model are the WMAP 
third year data together with the Hubble Space Telescope (HST) measure- 
ments (H — 72 ± 8km/s/Mpc [73]) and a top hat prior on the age of the 
universe between 10 Gyrs and 20 Gyrs. The resulting marginalised posterior 
probability distributions for the base cosmological parameters are represented 
in Fig. |TD] They are not significantly affected by the various prior choices on p 
and their corresponding mean values and confidence intervals are compatible 
with the current state or art [TJ. 

The constraints obtained on the small field inflation parameters are showed 
in Fig. [TJJ As expected, the allowed range for the power spectra amplitude 
V* is narrow. Concerning the above panels their interpretation require some 
precautions. Indeed, the probability that p takes small values depends on our 
theoretical prejudice on how big the field expectation value of the inflaton may 
be. If p is allowed to be much greater than m pl then all p are equiprobable. 
On the other hand, if np cannot take values bigger than 10 then small field 
inflation models with p ~ 2 are disfavoured. The p posterior shows on its 
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Fig. 10. Posterior probability distributions of the base cosmological parameters in 
the small field ylCDM inflation model. 



own that, independently of the p values, nfi > 10 is slightly preferred by 
the data. Since these posteriors are marginalised over the other parameters, 
the previous statements are necessarily robust with respect to any reheating 
model, in the framework of our modelisation. 

But more than being a nuisance parameter, Fig. II II shows that lni? is also 
mildly constrained by the data: the probability distribution of In i? has a lower 
bound slightly above the prior lni? > —46 given by Eq. (|65|) . Although this 
is not obvious, the upper cut-off seen in the lni? posterior comes from the 
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Fig. 11. Marginalised probability distributions for the small field inflation param- 
eters. The results coming from the prior k/x 6 [0.1, 100] are represented by solid 
black lines whereas, when they differ, the posteriors associated with the K[i prior in 
[0.1, 10] are plotted as dashed red lines [llj . 



upper bound of Eq. (|65j) (see Ref . [TT] for a more detailed discussion) . For K[i 
in [0.1, 100], we finally obtain at 95% confidence leveQ 

lni?>-34. (66) 

Plugging this inequality into Eq. (|64p constrains some properties of the 
reheating era. This result can be understood by looking at Fig. [5] Since vary- 
ing the reheating properties allow the observable window to move along the 
inflaton potential, it is not surprising that some part of the potential may be 
preferred from a data point of view. For the small field models, a more involved 
analysis would show that the spectral index of the scalar power spectrum gets 
away from what is allowed by the data when In R becomes too small [TT] . This 
is precisely why the current WMAP data lead to the bound ([6"6"| . 



7 The dependence of the hiR posterior distribution with respect to the Kfi prior 
disappears as soon as Kfj, is allowed to be greater than 10. 
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5 Conclusion 



As a conclusion, we would like to discuss the future directions associated with 
the possibility of constraining some basic properties of the reheating era with 
the CMB data. 

In the small field vlCDM inflation analysed in the previous section, the 
bound found in Eq. (|66p can be further explored by being more specific on 
the way the universe reheated. As a toy example, if one assumes that reheating 
proceeded with a constant equation of state P — wp, then Eq. (|64[) simplifies 
into 

1 — 3w , 4 s 1 + 3ui / 4 N , . 

ln R = 12TI^ ln {K Prch) + 6T6^ ln ( " Pcnd) ' (6?) 




Fig. 12. One and two-sigma confidence intervals (solid contours) of the two- 
dimensional marginalised posteriors (point density) in the plane [ln 7?, ln(fc 4 p en d)] 
for the small field models. The parameter k\i varies in [0.1, 100] from red to green. 
The four panels correspond to the situation in which the universe reheated with a 
constant equation of state P = wp. In each panel, the solid lines correspond to dif- 
ferent values of the reheating temperature (1/4) ln(/t 4 p re h) ranging from —45, —35, 
-25, -15 to -(1/4) ln( K 4 pond). 
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In the parameter plane [In ii, In (k p e nd)]j for a given value of the reheating 
energy, Eq. (|67p corresponds to a straight line. In Fig. five of these lines 
exploring the range p nuc < /3 ro h < p onc j have been superimposed to the two- 
dimensional probability distributions associated with the small field model, 
and this for four different equations of state having w > — 1/3, «; = 0, ui = 1/3 
and w < 1, respectively. 

At can be seen on the top left frame, for w > —1/3, low values of the 
reheating temperature are out of the confidence contours. From a more robust 
analysis using importance sampling, one would find in that case that, at 95% 
of confidence, p re h > 2 TeV [TTj . Of course this bound is not really impressive 
and close to the limits already set by BBN, moreover, it holds only for small 
field models and a quite extreme equation of state. However, is shows that it 
is already possible to get some information on the reheating era in a given 
model of inflation from the CMB data only. This is precisely on this point that 
particle physics models may be decisive. Indeed, once the inflaton couplings to 
the other particles are specified, the properties of the reheating era are fixed 
and a given particle physics model would appear as one curve parametrised 
by some coupling constants in the plane [lni?, ln(K 4 /j en( j)] of Fig. [T^l With the 
incoming flow of more accurate cosmological data, this may be an interesting 
way of constraining inflation as well as particle physics at very high energy. 
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